A STEADY STATE SOLUTION TO A MORTGAGE PRICING PROBLEM 
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Abstract. This paper considers a mortgage contract where the borrower pays a fixed mortgage 
rate and has the choice of making prepayment. Assume the market interest foUows the CIR 
model, the problem is formulated as a free boundary problem where the free boundary denotes 
the level of market interest rate at which it is optimal for the borrower to make prepayment. 
Here we focus on the infinite horizon problem. Using variational method, we obtain an analytical 
solution to the problem, where the free boundary is implicitly given by a transcendental algebraic 
equation. 



1. Formulation of the problem 

We consider a mortgage contract where the borrower pays a fixed rate of c {year~^) to the 
lender. In reality this mortgage rate is implicitly represented by a continuous payment of m 
%/year. At each time t when the contract is effect, the borrower has two choices: to continue 
the mortgage by paying mdt for the next dt period or to close the mortgage by paying off all 
the loan balance M{t), where the loan balance M{t) is determined by 

(1.1) ^ = -m + cM(t). 

When the contract duration T is given and M{T) = specified, the above ODE has a unique 
soultion 

(1.2) M(t) = -(l-e(-^(^-*))). 

c 

Here we assume the borrower always has sufficient amount of capital. The borrower chooses 
not to pay M{t) even though he is financially capable to do so if the expected future market 
return from an equal amount of investment is higher enough. On the other hand, if the expected 
future market return from an equal amount of investment is lower enough, he should choose to 
settle M{t). From the lender's point of view, the value of the contract V, as a function of time t 
and market interest rate x, is determined by the market interest rate x. The higher the market 
interest rate x is, the lower the contract value V, but V shall never be lower than 0. The lower 
the market interest rate x is, the higher the contract value V, but V shall never be exceed 
M(t) since the borrower has the choice to settle the loan once V reaches M{t). From standard 
mathematical finance theory, one can find the value of the contract V{x^ t) and the optimal level 
of market interest x = h{t) at which the borrower should make prepayment of M(t) by solving 
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the following free boundary problem: 



(1.3) 



' L{V) =m, 

c 

dV 

Vix,0)=0, 
[ h{0) = c 



for x>h{t),t>0 
for X <h{t),t > 

for all X > 



where t, for mathematical convenience, is defined to be the time to expiry of the contract, h(t) 
is the unknown free boudary to be determined together with V, and the differential operator L 
is defined as 



(1.4) 



L(F) 



dV 
'dt 



2"" dx^ 



k{e 



^9y 

x)— + xV 
ox 



The differential operator in the the system, referred as Kolmogorov equation, which can also 
be derived from Feymann-Kac Theorem [131 l3t [23]. Because of the important role played by the 
mortgage securities in real economy, there exists a considerable literature (see [H El [191 [3 [H] i 
for instance) dedicated to the topic, mostly of which have studied the problem from option- 
theoretical point with relatively less rigorous mathematical proofs. In recent development, a 
free boundary approach was introduced in [5l [10] to a similar problem. In particular D. Xie 
et al. (2007) have formulated the integral representations of the problems and proposed a 
Newton's iteration scheme under the assumption that the underlying interest rate follows the 
Vasicek model [5]. In this paper, we shall study the same type of mortgage contract with 
market interest rate following CIR model instead of Vasicek model. We use CIR model because 
it is observed that Vasicek model allows negative interest rate, which contradicts the empirical 
statistics from market [l]. In this paper, we focus on the infinite horizon (steady state) of the 
problem. 



2. The Infinite Horizon Problem 

Without loss of generality, we assume m = c, thus we derive the following infinite horizon 
problem of the original system. 



(2.1) 



„ ,2ke 2k ^ , 2 2c . * x 

XV + ( 5 n^W = X e [R , OO) 



a2 



V{R*) = 1 
V^{R*) = 
V{x = oo) = 



3 

where R* is to be determined together with V. To solve this free boundary problem, we let 

V{x) = e^^n(z) 



A 



2Vk^ + 2^2 

PX = 7i X 

k - + 2o-2 



The above infinite horizon problem is transformed into the following one. 

zu" + - z)u' - ^(1 - , ^ )u = -_^=e^^"i7^^', z e {z\oo) 
u{z*) = e'^~^^^^^'* 

u'(z*) = (- )r^2-,,/^2^,^2^^* 

^^^^ ' ^2 2VFT2^^ 

, ep^u(2) =0, z — > oo 

We first solve the ODE in z , disregarding any boundary conditions. The homogeneous equation, 
as a standard confluent hypergeometric equation, has two linearly independent solutions 

,^,k9. k . 2ke ^ 

17 V A;2 + 20-2 cr^ 
fc0 /c , 2ke 

U2 = U{ — {1 - ), — , Z) 

V A;2 + 20-2 cr^ 

one condition that both ^(1 — ^^2^2 2 ) ^'^'^ ^ ^^'^ strictly positive, which is certainly satisfied 
for this problem. And the Wronskian of these two linearly independent solutions is calculated 
to be 

r(^) 

W{M, U)iz) = 

By the standard variational method, we can find one particular solution to the inhomogeneous 
equation, and together with above two linearly independent solutions, we obtain the general 
solution of the ODE for u 

u{z) = c\ M{a, 7, z) + C2U (q, 7, z) 

r Mia,^,OU{a,^,z)-U{a,^,OMia,^,z) c ^^-(l+ ^y.^V^ 

Jz -M^-^ Vfc2 + 2a2^ 

and correspondingly, the general solution to the ODE in (j2.ip is given by 

V{x) = e^^{ci M{a, ■y,px) + C2U{a, 

(2-2) M{a,-f,OU{a,-f,px) - U{a,-f,OM{a,-f,px) r Ce~^^^ 
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where 



ke k 



7 



P 



2k9 

?2 



Next we need to decide, by the conditions at free boundary and infinity, the values of two 
constants in the general solution as well as the unknown infinite horizon z*. We firstly convince 
ourselves that ci = by investigation of the asymptotic behavior of the M(a, 7, z) and U (a, 7, z) 
for large z. Recall that 

M(a,7,z) ~ ^^z°-^e^ z ^ +00 

r(a) 

C/(q,7,2;) ~ x~", z ^ +00. 

If we let z — > 00 in the expression of u[z), the definite integral term vanishes, i.e., the U part 

contribution goes zero, but the M part go to e^. If we multiply them with , the magnitude 

contributed by the definite integral or the U part still go zero, which is the desired, because 

-2 

A < 0. But now the contribution of M will be multiplying a nontrivial polynomial in 2;, 

which clearly goes to infinity. This implies that M(a, 7, z) should not be included in the general 
solution. Now we only have two unknowns ci and z* . It is easy to use the fact that 

u{z*) = 2VPT5^^ 

u\z*) = (- ^ )e^^'2Vw)^* 

to get the implicit solution of z* : 

U{a+,'y + l,z*) 

^^^"JT^fep^^* _L c r°° M {a,'y,^)U (a,7,2*)-[/(a,7,g)M(a,7,2:*) 

g-7e 2\/fc^ + 2CT2 

^^^~ 2Vfc2+2^2 )^* , c r2(a) TOO A4'(a,7,g)t/(a+l,7+1.2*)-C/(«,7.0M(a+l,7+l,2*) ^c 

g-7e 2\/fc2+2CT2 

And the x* in original system can be simply recovered by x* = pz* . 
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